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Abstract 

Under the quenched approximation, we perform a lattice calculation for the 
mass of the ground 4 ++ glueball state in E ++ channel on a D = 3 + 1 
lattice. Our calculation shows that the mass of this state is Mc(4 ++ ) = 
3.65(6) (18)GeV, which rules out the 4 ++ or mainly 4 ++ glueball state inter- 
pretation for £(2230). 



1 Introduction 



As the experiments reported there maybe exist one candidate for glueball state, 
namely, £(2230), which has been seen mainly in J/ty radiative decay with a number 
of decay channels. But its spin has not been determined experimentally: we don't 
know whether its spin is 2 ++ or 4 ++ . Therefore, one hopes to calculate the masses 
of 2 ++ and 4 ++ glueball states on lattice for the spin interpretation of £(2230) since 
lattice QCD is regarded as the most reliable theory to study hadron spectra. 

Meanwhile, lattice simulation has achieved tremendous success in the study of 
glueball states. For instance, some authors show mass spectral! for glueball 
states. But, so far we do not know how to calculate the mass of 4 ++ glueball state 
assuredly on D = 3 + 1 lattice. On a D = 2 + 1 cubic lattice, in Ref. ||, Johnson 
and Teper performed a calculation for these states. Nevertheless, there are some 
disputes about the interpretation of the calculation, i.e., one does not know whether 
the mass calculated in Ref. || is really the mass of 4 ++ glueball states. In addition, 
as authors consider, one can hardly extent the construction of the operator in Ref. 
U into the D = 3 + 1 lattice. 

Therefore, from the point of views of both experiment and theory, it is important 
to calculate the 4 ++ glueball mass on lattice. Basing on the connection between 
the continuum limits of the asymptotic expansion of the chosen operator and the 
construction of the operator with definite J PC , we have developed a new procedure^ 
which can calculate masses of glueball states with arbitrary J PC . According to this 
procedure, a lattice calculation for the mass of the ground 4 ++ glueball states is 
performed in this paper by the construction of operators, corresponding states of 
which are in E ++ channel in the simulation. The leading term in expansion of 
this operator belongs to 4 ++ representation of SO(3) PC group. Therefore, from the 
asymptotic expansion of this constructed operators, we assure that, in continuum 
case, the spin of corresponding states is 4 ++ . Then, the mass of 4 ++ glueball state 
was obtained. We found that the calculated mass is prominently different with that 
of 2 ++ glueball states and it is about 1.5 and 2 times of that of tensor and scalar 
glueball state respectively. We should argued that our results rule out the possible 
interpretation of the ground 4 ++ glueball states for £(2230). 

This paper is organized as follows. In section 2 we show how to construct 4 ++ 
glueball operators in E ++ channel. We should also present our results and conclu- 
sions in section 3 and section 4. 



2 The Construction of the Operator Correspond- 
ing to 4 ++ Glueball state 

As we know, a definite state \ip > is generated by a current o acting on vacuum 
|0 >: 

!</>>= o|o>, (i) 



and the character of \ip > can be described by the current o. Therefore, as mentioned 
in 0, to calculate the mass of glueball state with spin J PC = 4 ++ , we should carry 
out three steps. First, we should write out nonet currents which transform as the 
representation 4 ++ under the SO(3) PC group and decompose them into irreducible 
representations , E ++ , 7\ and T 2 ++ of the group PC ( the cubic group with 
parity and charge conjugation transform ) in the subduced representation 1 . Then 
we construct corresponding operators on lattice. Its continuum limits should be just 
those currents mentioned above. Next we should also use the variantional procedure 
to extract masses of the ground 4 ++ states and excited 4 ++ states. 

We shall calculate the mass of the ground 4 ++ glueball states in this paper. 
Among nonet currents, there is only one type of such nonet current, the mass di- 
mension of which is the lowest one. This nonet current is of the form Tr(BDDB) 
and we expect that it should give the most contribution to the ground states with 
spin J PC = 4 ++ due to dimensional analysis. Therefore, we will only consider such 
nonet current and the construction of corresponding operators on lattice in this pa- 
per. Apparently such disposal will bring up some errors and its affectation will be 
considered in the error estimate in section 3. 



The Nonet Current and Its Decomposition To simplify, we denote magnetic 
field B x{yt z) = z) ^- by £ 1(2: 3 ) and covariant derivative D x ^ z) by D 1(2 , 3 ) 

(Di- = di — i[Ai, •]) respectively here. Obviously one can treat B^s as bases of the 
spin J = 1, and B± = T^(Bi ± iB 2 ) and B = B 3 as three standard basesQ. 
It is similar for covariant derivatives ZVs, which can also be combined into three 
standard bases: D± = T^(Di ± iD 2 ) and D = D 3 . 

Since spins of magnetic field B and covariant derivative operator D are both 
1, the nonzero gauge invariant nonet current of the mass dimension 6 is sole, so 
that according to C-G coefficients, one can write these 9 gauge invariant currents 
unambiguously in the standard form up to full covariant derivatives (The following 
Cj's(z = 4, • ■ ■ , —4) are eigen-currents of operator J 3 , i.e., J?,(ci) = i Cj.)0: 

c 4 = Tr [£+£>+£>+£+ + DG\, (2) 
c 3 = Tr[B + D + D B + + B + D + D + B + DG], (3) 

c 2 = -^Tr[(B + D + D_B + + B + D D B + + B + D + D + B_ 

+B D + D + B + 2B + D + D B + 2B + D D + B )}, (4) 

ci = -^Tr[£ + L> £>-£ + + B D + D + B_ + B + D + D„B 



1 The subduced representation used here is the representation obtained by trivially embedding 
the O group into the SO(3) p group. 

2 Here B± and B is the eigen-current of the operator J 3 , i.e., J^B^) = ±B± and J^Bq) = 0. 

3 Phenomenologically, we can regard this nonet current as 5 D^ multi-current and the corre- 
sponding state as 5 Z?4 multi-state. 



B + D_D+B + 2B + D D B + B + D + D B_ 

+B + D D + B_ + 2B D + D B + DG], (5) 

c = — L=Tr[fi + £L£L£ + + B-D+D+B-] 
2 Tr[B + D + D_B_ + B + D_D + B_] 
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+ -J=Tr[5+£L A,£o + B + D D_B 
v 70 



+B_D + D B + B D + D B_ 

+B + D D B_ + B D + D_B + B D D B + DG], (6) 

c_i = -^ r TY[B_D D + B_ + B D_D_B + + B_D_D + B 

B-D + D_B + 2B^D D B + B_D_D B + 
+B_D D_B + + 2B D_D B + DG], (7) 
1 

c_ 2 = -^TV[( J B_Z)_D + 5_ + B_D D B- + B_D_D_B + 

+B D_D_B + 2B_D_D B + 2B_D D_B ) + DG], (8) 
c_ 3 = Tr[fi_D_D B_ + B_D_.D_.Bo + DG], (9) 
c_ 4 = Tr[B_D_D_B_ + DG]. (10) 

Here is some full covariant derivative terms. In fact, one can get c_j(i = 
1, 2, 3, 4) easily by the swap + <-> — in Cj. 

Then we turn to the PC group, the finite subgroup of SO(3) PC group, to 
form the bases of subduced representation with nonet current mentioned above. As 
shown in Ref. 0, the subduced representation is reducible and one can reduce 
this representation into irreducible ones Af + , E ++ , T x ++ and T 2 ++ , the dimension 
of which are 1, 2, 3, 3 respectively. For our aim we only show the result of the 
decomposition into E ++ here. 

According to the procedure introduced in Ref. QlO | , we can get two bases of the 
representation E ++ 

ei = -^(c 2 + c_ 2 ) 

oc Tr[Bi{D\ - Dl)B l - B 2 (D 2 2 - D 2 3 )B 2 + B 3 (D 2 2 - Dj)B 3 
+AB 2 D 3 D 2 B 3 - 4B 1 D 3 D 1 B 3 + DG]; 

62 = ^ ( \/l c4+ /l c - 4 -/! co) 

oc TrfB^D 2 + D 2 - 2D 2 2 )B 1 + B 2 {D\ + D\ - 2D\)B 2 
+B 3 (D 2 + D\ - 2D\)B 3 + AB\D 3 D\B 3 

+AB 2 D 3 D 2 B 3 - %B 1 D 2 D l B 2 + DG). (11) 
These two currents are combinations of following three linear dependent ones: 



e« = Tv[B 2 {Dl-Dl)B % -Bz{Dl- 


-D 2 1 )B 3 + B 1 (D 2 3 


-Dl 




+4fi 3J D lJ D 3 5 1 - AB 2 D l D 2 B l 


+ DG); 






e< 2 ) = Tr[B 3 (Dj-D 2 2 )B 3 -B 1 (Dl 


-D 2 2 )B 1 + B 2 (Df 


-Dl 


B 2 


+AB 1 D 2 D 1 B 2 - AB 3 D 2 D 3 B 2 


+ DG}; 






= TiiB^Dl-DDB.-B^Dl- 


- Dl)B 2 + B 3 (D 2 2 


-Dl 


B 3 


+4B 2 D 3 D 2 B 3 - AB X D 3 D X B 3 


+ DG). 







(12) 



These three currents are of equivalence up to a cyclic permutation 1 — > 2 — ► 3 — ► 
1 of superscripts and subscripts in Eq. (|T2|). 

Then we consider how to construct operators on lattice, which are in the repre- 
sentation E ++ and their continuum limits are just eW(i = 1, 2, 3). 



The Construction of Operators on Lattice In this paragraph, we assign that 
(i, j, k) is even permutation of (1, 2, 3), i.e., (i, j, k) = (1, 2, 3), (2, 3, 1) or 

Using the link variable as defined in the discussion of improved actional H : 



U(n,i)=Texp( dtA^an + it) 



(13) 



we denote Wilson operators by figures Q. For instance, we denote the plaquette 
operator 



n) 



ReY,n[U(n,i)U(n+lj)U-\n + j,i)U-\n,j)] (14) 



by the figure XX J | ji ) for simplification. Here the summation is only over 

an I — £. 
all the spatial links in the same time slice. 

Let's consider the following five types of operators — d^: 

1. The first type of operators are d^\k = 1, 2, 3). They are from the following 
operators: 

frfc = e ijk + e jiki (15) 

£( i — 7 I + )• (is) 



where 



e (1) 



an 



t4 



an 



4 We use a notation, where j4j is an anti-hermitian, traceless N x iV-matrix. 



The expansion of the operator according to lattice spacing a is 



]T Tr[4 + 5a 4 (5| + B}) + ^(XID) + 5D 2 k + 12^)5., 



+ -^(17A 2 + 5Dg + 12^)5, - 2a*B i D i D i B 3 
+DG]+0(a 8 ). 
Here ^ = -\ £ e ij7£ F jfc and = diAj - djAi + [A+, Aj\. 

jk 

Then the operator d!jp is: 

4" = iP-ti" 

= £ TV[5o 4 (B? - B, 2 ) + ^B,(17DI + 5D] + 12Df)B 2 

n iZ 

-^(nDl + 5D* + YIB-DB, + a 6 B k (D* - Df)B k 
—2a 6 (BjD k DjB k - B^D.B^ + DG] + 0(a 8 ). 

2. Next, we consider operators d k 2 \k = 1, 2, 3). They are from operators 



(17) 



, (2) = (2) (2) 
k ^ijk ' ^jiki 



where 



c (2) 
-ijk 



E( 



an 



+ 



an 



(18) 

(19) 
(20) 



The definition and the expansion according to lattice spacing a of operators 
dtp are 

4 2) = &i 2) -&f 

= £ Tr[7a 4 (^ 2 - B*) + ^(WDj + + 24D*)B t 



-—Bj(19D 2 + 19D 2 k + 241%) + a?(BjD k DjB k - B t D k D t B k ) 



+DG]+0(a 8 ). 
Then we study operators df\ which are from operators 



(21) 



4 3) 



= e (3) 

ijk 

= E( 



an 



+ 




+ 



an 



/an ft 



(22) 



/an 



For d k we get 



d 



(3) 



E Tr [2a 4 (5j - 2?) + £^(0? + D\)B 3 - ^B % {D) + D^B, 

n 

,6 



+a (BjD k D 2 B k - B t D k D t B k ) + DG] + 0(a* 



Operators c?[ are from 



i>6 



(4) 



E< 



And again 



u 



(1) 



( 4 ) _ 

6 



= E Tr[^-(A 2 + D\)B 3 - B^D) + D\)B t 

n * 

+ DG]+0{a 8 ). 

5. At last we study the operator d k 5 ^ which are constructed by 



6f 



-ijk 



L. 



The construction is 



(5) 



fe (5) _ 6 (5) 



E Tr [2a 4 (E 2 - S|) + -^(Dj + 



(23) 



(24) 



(25) 



(26) 



--B J {D* + Dl)B J + DG] + 0{t 



(27) 



Therefore we let 

7(1) 



n n A n A 

= a 6 £ Trfi^A 2 - Dl)B { - B 3 \D) - D 2 k )B 3 + B fe (L>| - D 2 )5 fc 

n 

+4BjD k DjB k - ABiD k DiB k + DG] + 0(a 8 ). (28) 

While comparing Eq. ( |2"8"D and Eq. (|T2"|), we find that operators e k {k = 1, 2, 3) 
are just our aimed operators. Noticing that these operators are all 8-link, so that we 
need not consider the "tadpole" renormalization due to the mean field theory 0. 
We should utilize operators e k to calculate the mass of 4 ++ glueball states in this 
paper. 

Results and some discussions are shown in the forthcoming section. 



/(3) 



/(4) 



7(5) 



3 Calculation Results 



Under the quench approximation, we perform our calculation on an anisotropic 
8 3 x 40 lattice with improved action introduced in Ref. [O]: 



S 



11 



FiO 



+ 



3u 2 u? 



sir 



i2 U y t 



h 



(29) 



where (3 = Q/g 2 , g is the QCD couple constant, u s and u t are mean link renor- 
malization parameters(we set u t — 1), £ — a s /a t = 5.0 is the aspect ratio, and 
Q sp includes the summation over all spatial plaquettes on the lattice, Q tp indicates 
the temporal plaquettes, Q sr denotes the planar 2x1 spatial rectangular loops and 
Q s tr refers to the short temporal rectangles( one temporal and two spatial links). 
More detail is given in Ref. |12| . In each f3 calculation, we set 2880 sweeps to make 
configurations reach to equilibrium and then make total 8100 measurements in 90 
bins ( We perform measurement one time every four sweeps). The other calculation 
parameters, such as u A s) r s /r and a s , are as the same as those in Ref. H. 
Our results in different j3 are shown in Table 1: 



p 


1.7 


1.9 


2.2 


2.4 


2.5 


4 ++ 


1.61(1) 


1.35(7) 


0.97(4) 


0.79(3) 


0.70(5) 


2++ 


1.019(3) 


0.95(1) 


0.71(2) 


0.548(6) 


0.519(4) 



Table. 1 4 ++ glueball energy moat for each (3. The numerals in the brackets are 
error estimates( The data for 2 ++ glueball is from Ref. @.). 

Now we discuss a little about the error estimate. 

As discussed in many papers, one may regard that the chosen improved action 
breaks the rotation symmetry up to 0(a 4 ) in error estimate. But since the mass 
dimension of terms which break the rotation symmetry up to 0(a A ) is eight and the 
lowest mass dimension of our chosen operators in Eq. ( |2"%| ) is six, the up limit of 
precision here is 0(a 2 ). 

On the other hand, among the currents with mass dimension 7 or 9, there is no 
current with P = +. Therefore, the precision calculated here is actually 0(a 2 ) and 
the systematic error can be written out as the form C2a 2 + c±a\ + • • •, which should 
be used to fit our simulation data. Simulation results and fitting curve are shown in 
Fig. 1. The bottom curve and data@ in Fig. 1 are for the 2 ++ glueball state which 
also belong to E ++ channel. 

From the data and fitting curve, the statistical error of the mass measurement is 
0.044Ge\A As argued in ref. Q, there are 1% systematic error dur to aspect ratio. 
In addition, our approach leads into about 0.5% systematic error. Therefore, the 
total systematic error is about 1.1%( 0.041Gel^). So the mass of 4 ++ glueball states 
is 3.65(6)Gel / . Including the uncertainty in r^" 1 = 410(20)Mev, our final results 
are: M(j(4 ++ ) = 3.65(6) (18)GeV r which is about 1.5 times and twice of masses of 
the tensor and scalar glueball states respectively. 



4 Conclusion 



Through the construction of operators in E ++ channel which are complicated to 
some degree, we perform a calculation of mass of the ground 4 ++ glueball state 
under the quenched approximation. Due to the expansion of the chosen operator, 
we confirm that, in continuum, the mass calculated here is the mass of the ground 
4++ glueball state and it is M G (4++) = 3.65(6) (18)GeV. 

Apparently, our result rules out the J PC = 4 ++ glueball interpretation for 
£(2230) even by taking into account the opinion that the mass will be shifted by as 
much as 20% in going to full QCD®. While noticing that the common calculated 
mass for tensor glueball state is about 2AGeV under the quench approximation, we 
claim that our calculation supports the 2 ++ glueball interpretation for £(2230) if 
this state exists and it is confirmed glueball. 
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Figure caption 



Figure 1 The mass of 4 ++ glueball against the lattice spacing (a s /ro) 2 . Our 
fitting curve for the data of 4 ++ glueball state is m(4 ++ , a s ) = 3.65 — 1.22(a s /r ) 2 + 
2.74(a s /ro) 4 (wmt : GeV). In continuum case, the mass of 4 ++ glueball state is 
3. 65(4) GeV if we only consider the statistical error. The bottom data and curve: 
m(2 ++ , a.) = 2.417 + 0.783(a s /r ) 2 - 0.787 (a s /r )\unit : GeV) are for the mass of 
2 ++ glueball state which is shown in Ref. 0. 



